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Abstract 

We consider the electrostatic interaction between two rigid membranes, with 
different surface charge densities of opposite sign, across an aqueous solution 
without added salt. Exact solutions to the nonlinear Poisson-Boltzmann equa- 
tion are obtained and their physical meaning discussed. We also calculate the 
electrostatic contribution to the free energy and discuss the renormalization 
of the area per head group of the charged lipids arising from the Coulomb 

interaction. 
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I. INTRODUCTION 



Electrostatics often plays an important role in determining the structure of macro- 
molecules in aqueous solution, e.g. polyelectrolytes, charged membranes, and charged col- 
loidal particles. For many physical systems, where charge densities on the surfaces are equal 
and of the same sign, the Poisson-Boltzmann equation provides quantitative descriptions of 
their electrostatic interactions fflj. However, many biological processes involve charge densi- 
ties on the surfaces that are not equal and sometimes even have opposite sign 01. Examples 
of this situation include protein association with DNA and membranes; the interaction be- 
tween cationic liposomes and negatively charged cell membranes [[§. It also has significant 
biotechnological importance for the study of DNA association with artificial cationic lipo- 
somes Q. In a recent experimental work on two oppositely charged membranes 0, it is 
found that the membrane, due to its fluidity, adjusts its area per head in response to the 
electrostatic interaction. Similar conclusion has been reached by Radler et. al. || on DNA- 
cationic liposome complexes. Motivated by these experiments, in the present paper we study 
a model of two oppositely charged membranes in electrostatic interaction and quantify how 
the area per head renormalizes by using the Schulman-Montagne condition j7| for membrane 
self-assembly. 

We consider a system composed of two oppositely charged rigid membranes || separated 
by a distance L across an aqueous solution with dielectric constant e (see Fig. [l]). One 
carries positively charged lipids with magnitude q per head group and area per head E+. 
The other membrane is composed of negatively charged lipids of charge — q with E_ area 
per head. Without loss of generality, we assume here that El 1 > E7 1 . We imagine that 
the two membranes, each with its own counterions which we assume to be monovalent, 
initially infinitely far apart are brought into the vicinity of each other, where the electrostatic 
interaction dominates other interactions. Since the counterions can gain entropy by escaping 
to infinity, we only consider the positive counterions that are required to neutralize the 
system ||. We also assume that there is no salt in the solution in to order to focus on the 
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fundamental effects of the electrostatic interaction (see Conclusion for a discussion of the 
effects of added salt.) The charged lipids in the membrane are modeled as a fluid with an 
effective surface tension 7, which typically has a value of ~ 0.04 ksT/A 2 for bilayers; thus in 
addition to the electrostatic energy, there is a surface free energy of the form f3F s /N = 7E 
for each membrane, where N is the number of lipids, (3 = 1/ksT, T is the temperature, and 
ks is the Boltzmann constant. Note that energies are expressed in units of /cgT. We are 
interested in the values of E?j_ which minimize the total free energy - surface and electrostatic 
energy - as a function of L and 7. 

To illustrate the effect of the electrostatic interaction on the structure of membranes, 
let us consider the simplest case in which E + = E_ = E; in this situation there are no 
counterions between the charged membranes. The electrostatic contribution to the free 
energy can be calculated using Gauss' law to give JK| 



HP* = N p££) , (1) 
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where Ib = ~ 7 A, is the Bjerrum length for an aqueous solution of dielectric constant 
e = 80 and iV is the total number of lipids in the membrane. Minimizing F s + F e \ with 
respect to E, we find E* = ( ZL ^— J . Hence, the optimal area per head group depends 
on the square root of the distance separating the two membranes, and for a separation of 
~ 10 A, we find E* ~ 75 A 2 . In contrast, the electrostatic free energy for an isolated charged 
membrane with its own counterions is given by (Tl| 



/?F e ^-2AlnE, (2) 

to within an additive constant. Therefore, E* = 2/7 ~ 50 A 2 . Thus, the optimal area per 
molecule with electrostatic interactions may deviate substantially from that of an isolated 
charged membrane. Note that while the order of magnitude is the same in both cases, the 
functional dependences on the parameters are quite distinct. 

In order to take the counterions into account for unequally charged membranes, let us 
suppose that the negatively charged membrane is located at x = and the positively charged 



one at x = L. One constraint is that of charge neutrality - the number of counterions must 
be equal to the difference between the number of molecules on the two membranes: 

f L 11 

/ n{x)dx = — , (3) 

Jo h + 

where n(x) is the counterion density which is, in the mean field approximation, related to 
the potential (p by the Boltzmann factor: 

n(x) =n e- f3q<Kx) , (4) 

where q is the unit of charge and the prefactor no is fixed by Eq. Combining Eq. (f|) 
with the Poisson equation from electrostatics 

- V 2 0(x) = ^n(x), (5) 



e 

we arrive at the Poisson-Boltzmann (PB) equation: 

-^=4^oe^, (6) 

where we have defined ip(x) = f3 q<p{x). 

Equation (|6|) encapsulates a mean field approach to the many-body problem. It assumes 
that the counterions are point-like and collectively generate an average potential <f>(x) which 
governs how the counterions are themselves distributed. Furthermore, the PB equation, 
which neglects correlations among counterions, is valid only for sufficiently high temperatures 
or low surface charge densities fllj|l2|l. Therefore, within the mean-field approximation, our 
task is to solve the PB equation subject to the boundary conditions: 



and 



dip 
dx 



dip 
dx 



+^ (7) 



= (8) 



In the next section, we present the solutions to Eq. (|6]) and discuss the equilibrium 
configurations of the counterions. In section III we analyze the electrostatic contribution to 
the free energy and pressure of the system. In section IV, the phase diagram of the system 
is presented and followed by a discussion of the equilibrium value of the area per head. 



II. NONLINEAR POISSON-BOLTZMANN SOLUTIONS 



Equation can be solved using the "energy" method of classical mechanics, where we 
obtain a useful constant of motion, defined by 
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E = l^-j -8nl B n e-^ x \ (9) 

This constant can be physically interpreted as being proportional to the difference between 
the electrostatic stress and thermal pressure of the counterions (~ n(x)kBT). There are 
three cases to consider: i) E > 0, ii) E = 0, iii) E < 0. 
For E > 0, we have for the normalized potential 



ip(x) = In 
and the counterion distribution 



8irl B n _._ 1 _ 2 ( VE i 



E 



sinh (x — x 



(10) 



n{x) = —7=i — • (11) 

87r/ s sinh 2 (^|^) 

We choose the normalized potential to be zero at x — 0. This determines the value of 
no = n(0). Using the boundary conditions, one can show that the counterion densities at 
the surface of the membranes are given by: 

. . 2itIb E 

n{0) -li-m- B (12) 

and 

2nl B E 

n{L) = ^~^ (13) 



where E satisfies 



{± yf _ 4^( S+ -S-)^ VEL 



(14) 



For E = 0, we have for the normalized potential 

^p(x) = In 27iIbtio(x — x') 2 (15) 



and the counterion distribution 



nix) 



2nl, 



(x — x 



l\-2 



(16) 



The boundary conditions determine the value of L: 



2tt/ 



B 



(17) 



The density at each surface can be obtained using Eqs. (|T2|) and ( pT3|) by setting E = 0. 
For = — i£_ < 0, we have for the normalized potential 



ip(x) = In 



E_ 



cos 



(x — x) 



(18) 



and counterion density 



n\x) 



E_ 9 n/EI(x - x') 
sec 



8nl 



B 



(19) 



The boundary conditions give 

47T/ij(S_ 



E_ 



E_ cot 



(4tt/ b ) 2 



(20) 



Eqs. ( |i~2"D and (|T3|) are still valid with the replacement i£ — > —E_. 

We observe that for the dilute counterion limit, E > 0, the equilibrium density distribu- 
tion is essentially exponential. However, as the counterion density increases, the collective 
effect of mutual repulsion of the counterions leads to a sec 2 dependence of the counterion 
distribution as shown in Fig. |2|. Note that even for small separations, the counterions density 
is not uniform, in contrast to the case of electric double layers of like charges. 



III. THE ELECTROSTATIC FREE ENERGY AND PRESSURE 

The electrostatic free energy per unit area, which is the sum of the electrostatic energy 
and the entropy of the counterions, can be written as 

Pfei = J o dx \j^l n ( x ) [^O 17 ) -ln(rio^o)] | , (21) 

6 



up to an additive constant, where vo is the volume per counterion. From Eq. (pT|) or using 



a known expression for the pressure [12 



f3P = n(0) - (W>)787rZi 



we obtain 



P 



-E. 



(22) 



Here we have the simple result that the pressure is proportional to — E. Hence the solutions 
with E > and E < describe membranes that attract and repel each other, respectively. 
For the solution of E = 0, the membranes exert no net force on each other. Fig. |3| shows 
how the pressure varies with the distance. Using Eq. (|17]), we can determine the equilibrium 
distance L*\ 

* £+ - £- 

T * 1 

~~ 2ttI b 

Therefore, the result in Eq. (|22|) leads to the following picture. For large separation, L > L*, 
the counterion concentration is dilute and the electrostatic attraction dominates. On the 
other hand, the counterions are dense when the separation is small, L < L*. Hence thermal 
pressure dominates. When L = L* the electrostatic and thermal pressure balance, leaving 
zero net pressure. Similar results have also been obtained from numerical solutions to the 



PB equation for other systems (T^j. We note that since -qjj- > as suggested in Fig. |3|, 
the system is in a stable equilibrium at L = L*. 

Upon explicit evaluation of Eq. (|2T| ) and multiplying by the area of the membrane, the 
electrostatic free energy becomes: 



pF el /N + = -y(z) 2 z - 2 (A" 1 - 1) In z + A' 1 In fl - (A z y(z)) 



In 



l-{zy{z)f 



(23) 



where y± = 



E±L 



= L , and A = jjt. The function y(z) is defined as 



2 ' 
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y(z) 



Note also that y(z) ^ satisfies 



y+ if < z < (1 - A)" 1 
^ = (1-A)- 1 
iy- z > (1 - A)" 1 . 



(24) 



23/(2) = tanh[?/(z)] + A z y(z) (1 — z y(z) tanh[y(z)]) 



(25) 



In deriving Eq. (^), we have made an implicit assumption that the areas of the two 
membranes are the same, namely, iV + £ + = iV_£_. The validity of this assumption is justi- 
fied, since the edge effect, which arises when the areas of the two membranes are different, 
adds to Eq. (|23| ) a correction term of order 0(A~ 2 ). Therefore, in the limit of large sur- 
face areas, this contribution is negligible. Note that this assumption implies that the ratio 



A 



is fixed. Physically, the asymmetric parameter A gives the ratio between the num- 



ber of particles in each membrane. Equations (23), (|24|), and ( |25"D are the final results from 
which the equilibrium properties are derived in the next section. 



IV. EQUILIBRIUM PROPERTIES 

In this section, we determine the equilibrium value of £+, using Eq. (|23|) for the electro- 
static contribution to the total free energy, which is given by 

(3F tot /N + = 2 7 £+ + f3F el [L, S + ; A]/JV + , (26) 

where 7 is the surface tension. Let us first consider the case in which the system is at 
equilibrium (P = 0). By setting y{z) = in Eq. (p3|) , it can be shown that the optimal area 
per head group is given by 

A- 1 - 1 
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and the equilibrium distance can be calculated using Eq. ( |TTD to give 



(27) 



L* = 7- 1 (1 , A A )2 . (28) 
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Note that we must restrict A to the range of < A < 1, since the case of A = 1 corresponds 
to zero counterion density where the membranes attract for all separations. For the case of 
A = 0.1, we find ~ 200A 2 and L* ~ 5A. Therefore, the counterions have yet a stronger 
effect on the stretching of the membrane compared to the single charged membrane case 
where we have shown £ ~ 50A 2 . According to the discussion of Section [TTJ, Equation Q28D 
represents a line which separates the repulsive region and the attractive region in the phase 
diagram, which is shown in Fig. |j. 

Next, we consider two membranes separated by a fixed distance. Equation (|23D can be 
regarded as a function of only z with the help of Eq. fl24|) . Thus, we can imagine solving y(z) 
in terms of z, and by substituting the result into Eq. fl23|), minimization of the free energy 
can be carried out explicitly. Unfortunately, this cannot be done analytically in general. 
However, our estimate of L* above indicates that for L < L*, which is of atomic size, other 
effects such as thermal fluctuations and van der Waals attraction that have not been taking 
into account, may become significant. Therefore, we focus on the large separation limit, 
L >> L*, where y(z) » 1. With this approximation, we obtain 

S* + ^ ^ {l - A + ((1 - A) 2 + 47r/ B AV) 1/2 } . (29) 

Here, we observe that decreases as A — > 1 monotonically as shown in Fig. [5[ This is 
a reflection of the fact that the counterions have a strong effect on the stretching of the 
membrane. The presence of the counterions enhances the electrostatic energy; therefore, 
lowering of the charge density is energetically favorable. Hence the stretching of the mem- 
brane. For A = 0.1 and L = 10A, we find ~ 250A 2 which is five times the value of an 
isolated charged membrane! 

For completeness, we have performed the minimization of the free energy numerically for 
physically relevant separations (L > 2A). The results for the cases A = 0.1,0.5 are shown 
in Fig. |^. We see that Eq. (^) is indeed not a bad approximation for large distances. 
Furthermore, we observe that increases monotonically with distance for all values of A, 
as also suggested by Eq. (^). Note that even for L < L*, bears the same qualitative 
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dependence on A. 



V. CONCLUSION 

In this paper, we have considered two asymmetrically charged membranes of opposite 
sign. The counterion distribution is obtained by solving the PB equation. The membranes 
can either attract, repel, or exert no force on each other, depending on the difference in the 
area per head group of the membranes. We have also made an attempt to understand the 
collective effects of the counterions on the structure of the membranes. Due to the electro- 
static interaction of the counterions and the membranes, the stretching of the membranes 
becomes energetically favorable. In the high counterion concentration limit (A << 1) the 
area per head can reach a value few times as large as an isolated charged membrane. We 
note here that since the volume of a lipid bilayer is fixed, due to the incompressibility of 
lipid molecules, the stretching of the membrane implies a decrease in its thickness. This 
effect may be observed in real experimental settings by e.g. x-ray scattering. 

Finally, we wish to comment on the relevance of the model studied in this paper to other 
systems. First, for the case of permeable membranes, where counterions may freely permeate 
behind the charged membranes, the electrostatic interaction can be shown to be attractive 
for all separations, since the thermal pressure exerted on the membrane by the counterions 
is much weaker. In this case, we have verified that this problem is almost identical to our 
model for the solution E > 0. Therefore, the renormalization of area per head group as 
given by Eq. (|29D should hold at least qualitatively. Secondly, we also consider the presence 
of a finite amount of salt in our system. This problem may be described at the PB level by 
an equation similar to Eq. @, 

— — t = k 2 sinh?/' + AttIbUo e~^, 
dx z 

where k is the inverse of the screening length. The above equation can still be solved in 
principle but involves rather complicated mathematical functions. This is a subject for 
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further study. However, it can be inferred that our model corresponds to the case where the 
separation between the membranes is small compared to the screening length, i.e. L < n^ 1 . 
Typically, for 20A for 0.05 M of salt, our model should provide a good picture for the 

binding of oppositely charged membranes for a distance of the order of few A. For k^ 1 > L, 
the electrostatic interaction is exponentially screened. Nevertheless, according to Ref. J2] 
even with a finite amount of salt, the linearized PB theory also predicts regions of repulsion 
and attraction between membranes, similar to what we have found here. Therefore, we 
believe our model captures the essential physics associated with the electrostatic interaction 
of oppositely charged membranes. 
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FIGURES 




FIG. 1. Geometry of the problem: two membranes of unequal opposite charge, separated by a 
distance L with counterions between them. 
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FIG. 2. The counterion distribution for the case E < given by Eq. (|i~9|). It is not uniform for 
small distance in contrast to the case of two electric double layers of equal charge density of the 
same sign. 
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FIG. 3. A plot of pressure vs. separation distance for the case of A = 0.1 and £+ = 130A 2 



from Eq. (pi). 
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FIG. 4. The phase diagram as determined by Eq. (|2g|): a = 4-kIbLj vs. A. The region above 
the line is the attractive region of the membranes and the region below repulsive. 
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FIG. 5. The optimal value of the area per headgroup as a function of the charge asymmetric 
parameter A for L = 10 A. 
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